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New statistical method is proposed to coherently combine Baryon Acoustic Oscillation statistics 
(BAO) and peculiar velocity measurements exploiting decomposed density-density and velocity- 
velocity spectra in real space from the observed redshift distortions in redshift space, 1) to achieve 
stronger dark energy constraints, o(w) = 0.06 and a(w a ) = 0.20, which are enhanced from BAO 
or velocity measurements alone, and 2) to cross-check consistency of dark energy constraints from 
two different approaches; BAO as geometrical measurements and peculiar velocity as large scale 
structure formation observables. In addition to those advantages, as power spectra decomposition 
procedure is free from uncertainty of galaxy bias, this simultaneous fitting is an optimal method to 
extract cosmological parameters without any pre-assumption about galaxy bias. 
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I. INTRODUCTION 

The evolution of large scale structure, as revealed in 
the clustering of galaxies observed in wide-deep redshift 
surveys, has been one of the key cosmological probes. 
Structure formation is driven to grow by a competition 
between gravitational attraction and the expansion of 
space-time. This enables us to test our model of gravity 
at cosmological scales as well as the expansion history 
of the Universe. Although galaxies are not an unbiased 
tracer of the large scale matter distribution predictable 
from linear gravitational theory, various statistical meth- 
ods using galaxy redshift survey have been developed to 
probe the nature of cosmic acceleration such as dark en- 
ergy or modified gravity. Acoustic peak structure im- 
printed on density field inhomogeniety provides mapping 
to the large scale power spectra in transverse and longi- 
tudinal directions. Before decoupling, the photons and 
baryons are tightly coupled to form plasma, and over- 
density with baryon acoustic oscillation feature (BAO) 
develops through the competition between outward pres- 
sure support of radiation and gravitational attraction. 
After decoupling, as baryonic matter is influenced only 
by gravitational force at large scale, the baryons and dark 
matter instabilities form a configuration including pre- 
dictable BAO pattern extending from the sound horizon 
measured by CMB at decoupling epoch. This appears to 
be a standard ruler tracing the relation between redshift 
and expansion rate [ll-flij. Additionally, as galaxies are 
expected to act nearly as test particles within the cos- 
mological matter flow, the motions of galaxies carry an 
imprint on the rate of growth of large-scale structure and 
allows us to constrain cosmological models [l5l - |2]| . De- 
spite the intensive studies of both methods constraining 
cosmological parameters, the coherent way to combine 
both has not yet been fully optimized. 



Both different statistics can be combined coherently 
exploiting decomposed density-density and velocity- 
velocity spectra in real space from the observed redshift 
distortions in redshift space. The Fisher matrix exposi- 
tion was developed [22[ , assuming strict functional forms 
for the power spectra or allowing them to float freely. 
As expected the constraints are tightest when theoreti- 
cal investigations can provide good priors for the form 
of distorted power spectra in redshift space, but even 
relatively conservative assumptions suggest that percent 
level of decomposition should be possible with future 
surveys The cosmological distance errors are 

achievable using BAO preserved on decomposed density- 
density spectra, which is degraded a few factors of or- 
der from BAO of full observed power spectra. Decre- 
ment is compensated by adding cosmological constraints 
from the decomposed velocity-velocity spectra utilizing 
the growth of large-scale structure as imprinted on dy- 
namics of galaxies observed in large redshift surveys. 
Tighter constraints are achieved by this combination as 
well as dark energy constraints can be cross-checked by 
two different alternative approaches; distance measures 
and large scale structure formation observables. 
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The detailed formalism is presented in the next sec- 
tion. The Fisher matrix analysis to decompose spectra is 
briefly reviewed, then BAO statistics is redefined with the 
decomposed density-density spectra. Coherent formal- 
ism to combine BAO and decomposed velocity-velocity 
spectra is introduced, and advantages are discussed to 
gain stronger constraints on dark energy models with fu- 
ture surveys. For illustration, a fiducial ACDM cosmol- 
ogy is assumed with £7 m = 0.24, h = 0.72, n = 0.96 
and A 2 S = 2.41 x 10 -9 (in good agreement with a variety 
of observations) when computing specific predictions for 
future surveys. 
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II. COHERENT APPROACH TO COMBINE 
BAO AND PECULIAR VELOCITY 

The formalism is briefly reviewed to forecast errors 
on galaxy-galalxy density spectra, P gg , and velocity- 
velocity spectra, Pee (© is divergence of velocity map 
in unit of aH , = V • v/aH) from the observed spec- 
tra, P, of redshift surveys. We explain in detail how 
to estimate error forecast of BAO from the decomposed 
density-density spectra as well as how to combine de- 
composed peculiar velocity measurement coherently. 

A. Decomposition of the observed power spectra in 
redshift space 




The observed two-point correlation function in redshift 
space is decomposed into spectra of density fluctuations 
and peculiar velocity fields in real space 22j. The power 
spectra in redshift space, P, are given by, 

P(fc, /i, z) = {P gg (k, z) + 2 M 2 r(fc) [P gg (k, z)P ee (k, z)] 1/2 
+ li 4 Pee(k,z)}G FoG {k,n,a z ) (1) 

where (i denotes the cosine of angle between orienta- 
tion of two point correlation and the line of sight. The 
cross-correlation coefficient r(fc) is defined as r(fc) = 
Pg& I \/ Pgg Pee- The density and velocity fields are 
highly correlated for k < O.l/tMpc -1 thus wc assume 
that the density and velocities are perfectly correlated, 
r(k) ~ 1. The density-velocity cross-spectrum becomes 
the geometric mean of the two auto-spectra to leave only 
two free functions, P gg and Pee- Uncertainty in the ob- 
served rcdshifts is modeled by a line-of-smearing of the 

,222 

structure using fitting function GfoG = e where 
a z denotes one-dimensional velocity dispersion (FoG: 
Finger-of-God effect). 

The accuracy of decomposition of P gg and Pee out 
of P is estimated using Fisher matrix analysis determin- 
ing the sensitivity of a particular measurement. Fisher 
matrix for this decomposition, P^j c , is written as, 

dcc f dP(k) V eS (P) dP(k) d 3 k 



dp a P(kf dp/3 2(2tt) 3 
gg ' r ee- 



where p a = (P^ c , Pqq)- The effective volume V c g (P) is 
given by, 
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FIG. 1. Fractional errors of power spectra are presented in both 
panels. Upper panel: solid and dash curves represent fractional 
errors of decomposed Pgg C and observed P. Lower panel: solid 
curve represents fractional errors of decomposed Pee- Here o~(P) 
denotes the error on power spectra P, and cr(P)/P denotes the 
fractional error of it. 



where Pouter and Pinner denote comoving distances of 
outer and inner shell of the given redshift bin respectively. 
The weight function wf g(&, m) is given by 



WFoG(fc,Ai) = cxp 



(G 



FoG 
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where GfoG is the finger-of-god suppression factor and 
cT t h is a threshold value indicating our confidence in the 
accuracy of the FoG model. We fix er t h = 0.1 indicat- 
ing trust of clearance of non-linear effect up to 10% con- 
tamination of the factor beyond Kaiser limit. Derivative 
terms in Eq. ([2]) are given by, 



d\a.P{ki, jU, zj) 1 

dP^(k t , Zj ) = P{k l ^, Z] ) 

dlnP(ki,n,Zj) _ n 2 



1 + fi 



2 \ Pee(ki,Zj) 

Pgg(ki, Zj) 



Pgg{ki,Zj) 2 
I 1 



Pee(ki,Zj) 



(6) 



V eS (P) 



nP + 1 



(3) 



where n denotes galaxy number density, here n = 5x 
10 _3 ( h' 1 Mpc) -3 . Comoving volume, Purvey, given by 
each redshift shell from z = to 2 with Az = 0.2 {f 8 k y = 
1/2) is written as, 



47T 

fsky "g - {D, 
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B. Variances of power spectra 

The detectability of galaxy-galaxy density spectra in 
redshift space, P, is given by, 



1 



a 2 p (ki, Zj ) 



1 rkL 



^survey i^j) k^dkdfl 



(7) 

■j, j -i Jk' nlin [P(fc,/^) + l/n] 2 8tt 2 

where 'j' denotes redshift bin from z = to 2, and '%' 
denotes k space bin running from k — 0.01 h Mpc -1 to 
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0.5 h Mpc~ , i.e. integration of k in each ki bin from 
lower to upper bounds denoted by k l min and k z max respec- 
tively. 

The detectability of galaxy-galaxy density spectra in 
real space, P ggi can be estimated from decomposition 
Fisher matrix F^ c . The gg-component of the inverse 

matrix of P„£ c gives the variance of P gg c as, 



<Jpd, 



■■(h,z j ) = ^/(F^) gg {h, Zj ). 



(8) 



In the upper panel of Fig. [TJ solid curve represents frac- 
tional variance of P^ c at the same redshift bin 0.8 < 
z < 1.0. The decomposition procedure to transfrom mea- 
sured galaxy-galaxy density spectra in redshift space into 
real space degrades detectability, but is not very signifi- 
cant about a factor of 1.8. 

Variance of decomposed velocity-velocity spectra is de- 
rived from O0-componcnt of the inverse matrix of F^ c , 




,) = ^(jrdec-i )ee (9) 



In the lower panel of Fig. [I] solid curve represents frac- 
tional variance of P^e at redshift bin of 0.8 < z < 
1.0. As cth is set to be 0.1, the detectability is signif- 
icantly weakened at k > 0.1 /iMpc -1 . It is observed 
V-shape of variance pointed around threshold scale of 
k = 0.1 /iMpc -1 . Due to the short range of detectability, 
probing BAO structure using Pee proves to be challenge, 
but, as Pee is not a biased tracer, the measured ampli- 
tude can be used to constrain cosmological parameters. 

Threshold limit representing cutting-off scale of k due 
to the unpredictable non-linear physics, a^, reduces de- 
tectability of spectra shown in Fig. [TJ but differently af- 

The introduced <7th calibrates 



fects on Pgg C and 



r>dcc 

^ee- 



Fisher matrix integrand at the combination of k/i not k. 
The decomposition of P gg is less affected, as it is most 
determined at fi — > limit, while Pee is better decom- 
posed at fi — > 1 limit. Poorer detectability of Pee at 



k > 0.1/iMpc" 
fx — > 1 limit. 



is due to increasing FoG calibration at 



C. Formalism of BAO Fisher matrix 

Cosmological constraints arc forecasted from a mea- 
sured position of BAO in the three dimensional power 
spectra following the formalism of [ll|. The baryonic 
part of the power spectrum can be modeled as 



Pb(fc) 



sin kso 
ks 



■ exp 



Silk 



(10) 



where sq is the sound horizon at drag epoch and fcsnk is 
the Silk damping scale. The Silk damping scale can be 
accurately fitted by fc Silk = 1.6(ft b /i 2 ) a52 (ft m /i 2 )°- 73 [l + 
(lOAQrnh 2 )- - 95 ^- 1 hMpc' 1 . 



FIG. 2. Fractional errors on Dy(z) from Pgjj c 
rcprcsncdtcd by solid curve. 



BAO arc 



As in [Tl| , Eq. ([TO")) is multiplied by additional Gaussian 
functions to account for the erasure of information due 
to nonlinear evolution. The final Pb is given by 



PbO) = 



sinfcs „i/2 , M 



(11) 



where P0.2 is galaxy power spectrum at k = 0.2 /iMpc - , 
Aq is a normalization factor, S n i models the loss of infor- 
mation due to nonlinear growth, and GsAoik) is defined 



as. 



G BAO (k) = e^/feuO 1 



(12) 



Eq. (|10[) is compared with the fitting formula of Seo 
and Eisenstein to estimate Aq. The fiducial cosmol- 
ogy, as chosen in this research, has an analytical esti- 
mate of Aq — 0.46. As in [llj the numerical values 
S n ] = Eoff(5 m , where g$ m is the growth function normal- 
ized as gs m (a) = a at early time, and So = 11.0 ft. Mpc -1 
for the cosmology with erg = 0.8. When the physical 
value of the sound horizon is known to high precision 
from CMB measurements the error on sq is equivalent to 
the error on D v = (D^/H) 1 / 3 . 

Centroid approximation is available for P gg c , as it is 
decomposed into spherically symmetric Fourier space. 
BAO Fisher matrix from P^f° is given by, 



p2 
MD.2 



■i a P dBc(ki,Zj)' 
1 99 J 



■G 



BAO 



(h) , (13) 



where Nk is the total number of k bin upto k = 
0.1 /iMpc -1 . The variance of Dy is simply given by 

a(D^) = l/J~F$ of P^ c . In Fig. 03 fractional errors 

of Dy calculated from P^ g cc only are presented. 
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D. Tight constraints on dark energy from Pee 

The observed power spectra depend not only on dark 
energy parameters but the entire matter content and 
primordial power spectrum. We do not assume these 
quantities to be known. Instead, it is assumed that 
CMB data are available to constrain them. The CMB 
power spectra included in our analyses are the Ci spectra 
of temperature-temperature, temperature-polarization 
and polarization-polarization. Cosmological parameter 
space is given by (w, w a , u>b, w m , 9s, As, n s , z reion ). 

With CMB temeprature fixed at T C mb = 2.726 K, the 
sound horizon at last scattering surface, Ds{zi ss ), is given 
in Mpc unit by ui m and wj, measured by CMB acoustic 
peak structure. With the measured angular size of sound 
horizon at last scattering surface, the angular diameter 
distance to the last scattering surface is also given in 
Mpc unit. The geometrical factor determination in Mpc 
unit is converted into the Fourier space dimension of k 
in Mpc -1 unit. As primordial spectra are tightly con- 
strained by CMB physics, the pivoting scale of primor- 
dial spectra k p is given in Mpc -1 unit. In addition to 
this, as the shape of spectra in terms of k (transfer func- 
tion X$) is tightly determined by parameters measured 
by CMB physics, we express T$ in terms of k in Mpc -1 
unit. 

In redshift survey, measured redshift of objects is con- 
verted using distance in h^ 1 Mpc unit in which mapping 
is less affected by most uncertain distance parameter of 
the Hubble constant. The velocity-velocity power spec- 
tra are given by in terms of k in /iMpc -1 unit (2lj . 
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Pee(k,a) 



8tt 2 k 
25 #?Q2 



A 2 
As 



kh 



Tl(kh)gl{a) 



Variation of dark energy parameter impact on Pee(&; 
becomes complicated when we consider the combined 
constraints from CMB and redshift survey in which k de- 
pendent quantities are measured in different units. First, 
multiplication factor becomes k/HfH 2 n which changes 
due to floating of VL m and h. Second, primordial tilt fac- 
tor changes as k and k p arc determined by different units. 
Third, transfer function factor T$ is unmodified in terms 
of kh with varying dark energy parameters, as decay due 
to silk damping is well determined by cosmological pa- 
rameters measured at last scattering surface with CMB 
physics. But it is altered due to using unit of h" 1 Mpc 
instead of Mpc for distance measure of redshift survey. 
Finally, growth factor changes. 

In Fig. El total change of P e e(k = 0.86 h Mpc -1 , a) is 
presented with varying w (upper panel) and w a (lower 
panel). If degeneracy is not considered between w and 
w a , tight constraints on dark energy paramters are ob- 
served. Fourier space is split into 10 bins from k = 
0.01ft,Mpc _1 to 0.5/iMpc" 1 . This figure is a repre- 
sentation by one of 10 k bins chosen at 0.07 < k < 
0.1/iMpc -1 . 
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FIG. 3. Solid curves represent fractional error of Pqq at k range of 
0.07 < k < 0.1 h Mpc -1 in both panels. Upper panel: dash curves 
represent percentage variation of Pqq with different w from w = 
-1, w = _o.95,-0.96,-0.97,-0.98,-0.99 from top to bottom. 
Lower panel: dotted curves represent percentage variation of Pee 
with different w a from w a = 0, w = 0.15,0.12,0.09,0.06,0.03 from 
top to bottom. Here, k = 0.86 Mpc -1 for both panels. 



E. 



Coherent combination of Pg° c and Pqq 



To calculate the expected parameter errors we make 
a first order Taylor expansion of the parameter depen- 
dence including CMB, BAO and decomposed Pee- The 
expected parameter errors are easily calculated from the 
inversed covariance matrix of the Fisher matrix using this 
'linear response' approximation. The linear response ap- 
proximation is improved and susceptibility to numerical 
error is reduced with a careful choice of the parameters 
including 6 S . F™ B denotes a Fisher matrix of CMB 
physics equivalent to combination of distance measure 
at last scattering surface and primordial parameter con- 
straints, thus exploiting unlcnscd CMB power spectra 
only. We use Planck experiment specs at three frequen- 
cies of 100, 143 and 217 GeV. Each channel has the reso- 
lution of 9.1, 7.2, 5.0 arcmin , the temperature sensitivity 
of 5.5, 6.0 and 13 f^K and the polarization sensitivity of 
oo, 11, 27 iiK respectively. We calculate Fj^ B using the 
method developed in our previous paper [24|. 

Constraint on distance measure Dy by BAO is con- 
verted into cosmological parameter constraints using the 
following Fisher matrix formalism, 



F, 



D ,™ 



2 P 

^2 ^ 



i 



z i= 



dXr, 



pde 



8x n 



(14) 



where indices m and n denote elements in cosmolog- 
ical parameter space. Dark energy constraints from 
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FIG. 4. We present the fractional error dependence of dark energy 
constraints from BAO in the upper panel of Fig. [T] outer dash 
contour represents dark energy constraint from the fractional error 
of Pgg C (solid curve in the upper panel of Fig. [TJ, and inner dash 
contour represents dark energy constraint from the fractional error 
of P gg (dash curve in the upper panel of Fig. [I}. 



Fmn + F%^ B arc shown to be the outer dash contour 
in Fig.|4]and the dash contour in Fig.[5j <r(w) — 0.15 and 
a(w a ) = 0.37, which is weaker than dark energy con- 
straint from BAO of full observed spectra, a(w) = 0.10 
and <j(w a ) = 0.30 (the inner dash contour in Fig.0}. Here 
error estimation a denotes 68% statistical significance. 

pdec 

While D v 39 constrains dark energy by distance mea- 
sure, dark energy is probed by signature on the large 
scale structure growth of peculiar velocity measurements 
Pee- Shown in the bottom panel of Fig. Q] for fractional 
errors of Pee, the effective range of k is too narrow to 
observe the baryonic acoustic peak structure. However 
unlike P^ c , it is an unbiased quantity that is predictable 
with a given set of cosmological parameters. Fisher ma- 
trix is given by, 



F ee = 

mn 



EE 



1 



dP@e(ki,Zj) 



k — l Zi 



dx r 



FIG. 5. Constraints on dark energy parameters are presented; dark 
energy constraint from Pgg C BAO (dash contour), Pee (dotted 
contour) and coherent combination of both (solid contour). 



In Fig. [51 dotted contour represents dark energy con- 
straint from Pqq- Tight contraints achieved using Pqq 
is worth being investigated. At last redshift bin of 
1.8 < z < 2.0, dark energy constraints can be read from 
Fig. [3] about a(w) ~ 0.035 or cr(w a ) ~ 0.09 which con- 
firms the results from Eq. [15] If fitting is extended to 
all redshift bins then o{w) = 0.015 or a(w a ) = 0.04 (de- 
rived from Eq. ll5[) . This still agrees with that is observed 
in Fig. [3j But simultaneous determination of w and w a 
becomes less tight as a(w) = 0.18 and (j{w a ) = 0.50, be- 
cause of degenerate Pee variation in terms of w and w a . 
The results shown in Fig. [5] are given after full marginal- 
ization with all other cosmological parameters. 

Then coherent combination of BAO and peculiar ve- 
locity measurements can be simply given by 



F, 



total 



F„ 



pee , pCMB 

win ' mn 



(15) 



Solid contour of Fig.[5]rcprcscnts dark energy constraints 
from F^\ a(w) = 0.06 and cr(w a ) = 0.20, which is 
tighter than any other case considered in this paper. 



At scales greater than k = 0.1/iMpc - , the uncer- 
tainty increases due to the unknown non-linear physics 
causing the FoG effect and the violation of assumption 
of perfect cross-correlation between density-density and 
velocity- velocity spectra [H[ . It was shown that the de- 
composed Pee starts to be biased as scale approaches 
to k = 0.1/iMpc -1 . Cut-off scale needs to be intro- 
duced at fc cut = 0.1/iMpc -1 in order to skim away 
contributions contaminated from non-linear physics at 
k > 0.1/iMpc" 1 . 



III. CONCLUSION 

The known galaxy bias has always been considered an 
obstacle utilizing redshift surveys for cosmological pur- 
pose. If just peak structure is considered as a standard 
ruler at the given targeted redshift, cosmological infor- 
mation is extractable with less dependence on bias uncer- 
tainty. However enriched information of large scale struc- 
ture is not fully exploited. Previously, decomposed con- 
tributions of density-density and velocity-velocity spec- 
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tra were studied in the observed redshift distortion using 
anisotropic feature aligned along the line of sight. In this 
paper, it is proposed to coherently combine both decom- 
posed spectra for constraining cosmological parameters 
(still in bias uncertainty independent way) to enhance 
dark energy constraints. 

What is implemented in this new method is to cross- 
check constraints on dark energy in two different ap- 
proaches; distance measure by BAO and large scale stuc- 
ture growth factors by peculiar velocity measurements. 
There are many different theoretical models proposed to 
generate cosmic acceleration. It is essential to test those 
theories in many different observational methods. When 
we decompose spectra into density-density and velocity- 
velocity pairs, we can test theoretical models using geo- 
metrical tools as well as using large scale structure forma- 
tion. In addition, dark energy constraints are enhanced 
from performing BAO alone using density-density spec- 
tra before decomposition. In this sense, our approach is 
optimizing way to exploit redshift surveys for the purpose 
of constraining dark energy. 

It is important to test consistency between geometrical 
measures and structural observable. As general relativ- 
ity can be modified at large scales. One of vital mission 
of dark energy probe is to distinguish two different class 
of models; modifying energy momentum components or 



modifying gravity. Coherent motion measurements pro- 
vide us with a way of tracing Newtonian force sourcing 
dynamics of gravitational instability. Our new statistical 
methods assist us testing modified gravity type models 
too. 

It has not yet been fully tested whether the assump- 
tion of preservation of baryon acoustic peak structure is 
true or not, although baryon acoustic structure is not ex- 
pected to be lost in the decomposition procedure. It is 
necessary to test this assumption using N-body simula- 
tion. For decomposition of velocity measurements, it was 
tested whether there was velocity bias or not. However at 
the same time, it was known that velocity- velocity spec- 
tra are most sensitive to be extracted without systematic 
uncertainty. In the future, higher resolution experiments 
will be launched, and it would be interesting whether 
velocity bias is not induced yet. 
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